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Abstract -We calculate the thermoelectric figure of merit, zT = S'^GT/{ki + Ke), for p-type Si 
nanowires with axial Ge heterostructures using a combination of first-principles density-functional 
theory, interatomic potentials, and Landauer-Buttiker transport theory. We consider nanowires with 
up to 8400 atoms and twelve Ge axial heterostructures along their length. We find that introducing 
heterostructures always reduces S^G, and that our calculated increases in zT are predominantly 
driven by associated decreases in kj. Of the systems considered. 111 nanowires with a regular 
distribution of Ge heterostructures have the highest figure-of- merit: zT = 3, an order of magnitude 
larger than the equivalent pristine nanowire. Even in the presence of realistic structural disorder, in 
the form of small variations in length of the heterostructures, zT remains several times larger than 
that of the pristine case, suggesting that a^xial heterostructuring is a promising route to high-zT 
thermoelectric nanowires. 



In recent years, silicon nanowires (SiNWs) have been proposed for use as chemical sen- 
sors [U , photovoltaics P] and thermoelectrics [3113] . Hicks and Dressclhaus [5] first identified 
that NWs could be used to improve the thermoelectric figure of merit, zT, over bulk and 
two-dimensional superlattice (2DSL) values. While the dramatic increases that were pre- 
dicted have not been yet realised, much progress has been made: Refs. report zT ^ 1 
for SiNWs, a 100-fold increase over the value for bulk Si. 

The measure of the performance of a thermoelectric material is given by its figure of merit 
zT = S'^GT/ {ki + Ke), where S, G and T are the Seebeck co-efficient, electronic conductance 
and average temperature of the two contacts, respectively, and ki and Kg are the lattice and 
electronic contributions to the thermal conductance, respectively. zT may, therefore, be 
increased by designing materials that have either higher thermoelectric power factor S'^G, 
or lower thermal conductivity. The relatively high zT seen in recent experiments on SiNWs 
has been attributed to both of these effects: an increase in S resulting from enhanced phonon 
drag [3], and a decrease in ki resulting from surface scattering of phonons 

By analogy with 2DSLs, which show a reduction in k; compared to its value both in 
the bulk and in the alloy limit [7], super latticed NWs have been proposed as a possible 
route toward high zT thermoelectrics [8], through both enhancement of the power factor 
S^G and reduction of ki, as compared to pristine NWs. Experimental evidence for Si-SiGe 



(""^ Email: a.mostof iOimperial .ac.uk 



p-1 



Matthew Shelley and Arash A. Mostofi 



<111> <211> 





Fig. 1: Top: Cross-sections of SiNWs. Labels indicate the crystal direction of the longitudinal 
(z) axis (pointing into the page). Diameters (left-to-right): 0.78 nm, 1.02 nm, 1.44 nm, 1.14 nm 
and 1.06 nm. Middle: A single Ge heterostructure in a SiNW (Ge and Si atoms in magenta and 
blue, respectively). Bottom: A multiple heterostructure nanowire (MHNW) with an arbitrary 
distribution of Ge heterostructures. 



NWs [S] supports the idea that superlatticing results in a reduction in k;, although direct 
experimental evidence for increased S^G is still missing. 

In this Letter, we use a combination of first-principles density-functional theory simula- 
tions and calculations with interatomic potentials in order to compute zT for axially het- 
erostructured Si-Ge NWs within the coherent transport regime using the Landauer-Buttiker 
approach [TOIIII]. More specifically, we calculate zT at 300 K for thin (< 2 nm diameter), 
p-type (110), 111 and (211) H-passivated SiNWs (Fig. [TJ top) containing: (a) single ax- 
ial Ge heterostructures with lengths ranging from 0.4 nm to 4.3 nm (Fig. [TJ middle); (b) 
multiple Ge heterostructures of uniform length, distributed along the length of the NW 
either randomly, periodically, or as a Fibonacci chain (Fig. [TJ bottom); and (c) multiple Ge 
heterostructures whose lengths approximately follow a Gaussian distribution. 

Our calculations demonstrate that: (i) the introduction of a single Ge heterostructure 
in a SiNW can lead to a 3.5-fold increase in zT as compared to the equivalent pristine 
SiNW; (ii) this can be further enhanced to a 7.4-fold increase by introducing multiple Ge 
heterostructures and by controlling their spatial distribution along the length of the SiNW; 
(iii) this observed enhancement in zT is almost entirely due to a reduction in thermal con- 
ductivity rather than an increase in the power factor S^G - indeed, we find that introducing 
Ge heterostructures in an SiNW always results in a decrease in S^G; and (iv) introducing 
further disorder, in the form of a variability in the length of the Ge heterostructures within 
a SiNW, results in a decrease in zT as compared to the case in which the heterostructures 
are all identical, highlighting the importance of precise atomic-scale control that may be 
required in order to fabricate high-zT NWs. 

Our method (described in detail elsewhere [T^j), in which accurate yet compact model 
Hamiltonians of large-scale systems are constructed from first-principles calculations, enables 
us to study transport through meso-scale systems with modest computational cost: our 
largest simulations consist of a conductor region of length 116 nm (8432 atoms) coupled 
to semi-infinite leads. Our procedure is largely automated, which has made it possible 
to perform high-throughput calculations and undertake a comprehensive study of a large 
structural parameter space. With little modification, our general approach may be easily 
used to calculate transport properties in other quasi one-dimensional systems. 

We note that analogous approaches have been used recently for calculating electronic 
and/or thermal transport in large-scale ID systems, e.g., in SiNWs [T3J[T3| and carbon 
nanotubes [T5HT7] . 
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Fig. 2: Maximum thermoelectric power factor S^Gmax at 300 K for single Ge heterostructures in 
(110) (black solid lines), 111 (red dashed lines) and (211) (blue dot-dashed lines) SiNWs. Three 
diameters 0.78 nm (triangles), 1.02 nm (circles) and 1.44 nm (squares) are shown for the (110) 
direction. Pristine SiNWs are shown plotted as zero heterostructure length and pristine Ge NWs 
values are shown on the right. 



Electronic transport properties. — Starting from a plane-wave density-functional 
theory (PW-DFT) calculatioiill, a unitary transformation is applied to the extended ground 
state eigenfunctions in order to obtain maximally-localized Wannier functions (MLWFs) [50] 
and, hence, the Hamiltonian matrix in the basis of MLWFs. Due to the localized nature 
of MLWFs in real-space, the Hamiltonian matrix can be spatially partitioned and used 
in so-called "lead-conductor-lead" Landauer-Buttiker [TUIIII] transport calculations, using 
standard Green function techniques [THI1TH53] . In our case, the semi- infinite "leads" are 
pristine (H-passivated) SiNWs and the central "conductor" region comprises of SiNW with 
some axial distribution of Ge heterostructures. An example of such a system exhibiting a 
single Ge heterostructure is shown in Fig. [T] (middle panel). Although we begin with PW- 
DFT calculations with periodic boundary conditions, we determine electronic transport 
properties under open boundary conditions. Once the electronic density of states of the 
conductor and the transmission function 7~(e) are calculated, one can write I251[5B] G = 
e^Lnip), S = Li{^)/eTLo{^i) and Kg = ^{L2(/i) - [Li{^i)]'^ / Lo{^i)}, where 

L.„M = 2£rW(.-rt™(^^)^., (1) 

and f{e,^) = l/{exp[(e — p,)/kBT] + 1} is the Fermi-Dirac function at chemical potential 
fi. In this Letter, we focus on hole transport, so that we can associate fi with a carrier 
concentration that is driven by p-dopin^. 

As full DFT structural relaxation of our large heterostructured NWs would have been 
computationally intractable, the atomic configurations used for the electronic calculations 
were built by piecing together unit cells of pristine Si and Ge NWs whose equilibrium lattice 
parameter^ were calculated separately with DFT. For the smallest NWs considered, the 
results obtained from our approach had only small quantitative and no qualitative difference 
as compared to those from the equivalent fully relaxed structure [T^ . 

Fig. [5] shows the maximuno thermoelectric power factor S^Gmax for a range of Ge 

^We use the Quantum-Espresso package 1181 . the local-density approximation for exchange and correla- 
tion, norm-conserving pseudopotentials, a 400 eV energy cut-off for the PW basis set, and F-point sampling 
of the Brillouin zone. 

^This doping is 'artificial' in the sense that we do not directly include dopant atoms in our calculation. 

^These were 3.775 A (3.910 A), 9.224 A (9.497 A), and 6.470 A (6.692 A), for our largest diameter Si 
(Ge) NWs in the (110), 111 and (211) growth directions, respectively. 

■'As can be seen from Eq. 1^, the electronic transport coefficients are functions of chemical potential ^. 
Throughout this work, S^Gmax and zTmaic are defined to be the maximum values of the power factor and 
the figure of merit, respectively, as a function of fi. The maximal power factors shown in Fig.[2]are obtained, 
in all cases, for values of fj, within 30 meV of the Si valence band edge. 
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heterostructure lengths in (110) (black solid lines), 111 (red dashed line) and (211) (blue 
dot-dashed line) SiNWs. We note first that in no case does the introduction of an axial 
heterostructure result in an increase of the value of which is, at best, approximately 

the same as that of a pristine SiNW. The similarity of the results between pristine SiNWs 
and GeNWs is also interesting to note. Since all the NWs investigated only have a single 
channel that is available for conduction at the top of the valence manifold, wc confirm that 
the most important factor for S'^Gmax in quasi-onc-dimcnsional systems is the number of 
conducting channels at this edge [57] ■ 

The oscillations in Fig. [5] can be explained by a model in which the heterostructure 
is considered as a ID quantum potential well of width L, corresponding to the length of 
the heterostructure, and depth Vq, corresponding to the band offset between Si and Ge. 
The reflection amplitude of a wave incident on such a well vanishes when the well-known 
Fabry-Perot resonance condition is satisfied, qL = rnr, where q is the wavevector inside 
the well and n is an integer. For holes entering the heterostructure, this condition gives 
En{L) = —2^X2 ~^ '^here m* is the effective mass of holes and En{L) are energies at 
which resonances occur in the transmission. For the 111 NWs, which show the strongest 
oscillations in S'^Gmax, plotting the resonance energies against 1/L^ produces an excellent 
linear fit (not shown), with correlation coefficient = 0.998, giving m* = 0.28 me and 
Vq = 0.32 eV. Models such as this may be a useful additional tool for the optimization of 
heterostructure lengths in thermoelectric devices [28] . 

Phononic transport properties. — We determine the lattice thermal conductance 
Ki in an analogous way to the electronic conductance, by segmenting the system into a 
lead-conductor-lead geometry. Instead of finding the Green function of a Schrodinger-type 
eigenvalue problem, we determine the Green function that solves the eigenvalue problem 
relating nuclear displacements u to the dynamical matrix K and phonon frequency w: Ku = 
ijpu. Applying the thermal equivalent of the Caroli formula [21][29], we obtain the phonon 
transmission function, T{uj), in the limit of non-interacting phonons and coherent phonon 
transport. One can then write [30] 



Determination of the dynamical matrix using first-principles methods is computationally 
intractable for the large NW supercells with heterostructures that are considered here. To 
structurally relax the NWs and obtain their dynamical matrices, therefore, we use TersofF 
potentials [31] , which have been shown to give accurate values for lattice thermal conductiv- 
ities for thin pristine SiNWs, as compared to DFT calculations [35]. It is worth noting that 
the approach outlined above neglects Umklapp scattering, which would further decrease k/ 
at the temperature with which we are concerned (300 K). 

Results for ki arc shown in Fig. [3] (top panel) for the same single Ge heterostructure 
SiNWs discussed earlier. We find that, by introducing a single Ge heterostructure, the 
lattice thermal conductivity can be reduced by a factor of five for the 111 growth direction, 
as compared to the corresponding pristine SiNW, giving a value of 0.1 nWK^^. Reductions 
in the (110) and (211) direction are also significant (approximately a factor of four). We 
tentatively suggest that the longer unit cell in the 111 direction may account for the greater 
reduction seen in that growth direction. It may also be seen that, in the (110) direction, 
Ki increases with diameter as more phonon modes become available. We also observe this 
trend in larger diameter pristine SiNWs in the 111 and (211) directions (results not shown). 

zT for SiNWs with a single Ge heterostructure. — Fig. [3] (bottom panel) com- 
bines our results for electronic and phononic transport coefficients for the single Ge het- 
erostructure systems discussed above and shows our calculated values of zT^ax as a function 
of heterostructure length. It can be seen that heterostructured SiNWs in the 111 direction 
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Fig. 3: Dependance of m (top panel) and zTmax (bottom panel) at 300 K as a function of Ge 
heterostmcture length for (110), 111 and (211) SiNWs. Labelling is equivalent to Fig. [2] again 
plotting pristine SiNWs as zero heterostmcture length and pristine Ge NWs also shown on the 
right. 



display the greatest values of the figure of merit zT^ax — 1-4. Such high values, however, 
are not found consistently across the range of heterostructures studied. Such variations 
may limit the values of zT observed in realistic SiNWs since experimental control over het- 
erostmcture length is, currently at least, limited to length-scales comparable to, or greater 
than, the differences in length that are investigated here [Sj- In the (110) and (211) direc- 
tions, we find zTniax < 1, mainly due to the higher lattice thermal conductivities found in 
these systems. We note that, across the range of systems studied, the ratio of lattice and 
electronic thermal conductances, kz/ke, lies between 3 and 10, therefore, is the dominant 
contribution to the denominator of zT. This emphasizes the importance of reducing the 
lattice thermal conductivity for high zT NWs. 

zT for SiNWs with multiple Ge heterostructures. — Next, we consider much 
longer SiNWs with many Ge heterostructures along their length. Such multiple heterostmc- 
ture nanowires (MHNWs) are shown schematically in Fig. 2] These systems are too large 
for brute-force PW-DFT calculations. Instead, we use the Hamiltonian matrices of single 
heterostmcture calculations as "building-blocks" for constructing model Hamiltonians of 
much larger (up to ~ 8400 atom) MHNWs, with negligible loss of accuracy. Our approach, 
which relies on exploiting the nearsightedness of electronic structure that becomes manifest 
when Hamiltonian matrices are represented in a basis of MLWFs, is described in detail in 
Ref. [ig. 

Once the model Hamiltonian for the MHNW is constructed, the electronic transport 
properties under open boundary conditions are calculated in exactly the same way as de- 
scribed above for SiNWs with single Ge heterostructures. For the lattice thermal conduc- 
tivity, an analogous "building-block" scheme is used in which the (short-ranged) dynamical 
matrices of the single heterostructures that comprise the MHNWs are combined to construct 
dynamical matrices for the MHNWs. Under the assumption that phonons remain phase co- 
herent across the length of the MHNW, this dynamical matrix is then used to calculate the 
phonon transmission function T{uj) and, hence, the coherent lattice thermal conductance 
K™^ according to Eq. 

It is unclear whether the phase coherence length of phonons is comparable to the lengths 
of the MHNW systems that we consider (up to ~ 116 nm). Therefore, we also calculate 
the lattice thermal conductance in an ohmic regime, k°^™, in which the total resistance of a 
given MHNW is the sum of the thermal resistances of each individual heterostmcture that 
constitutes the MHNW [HI]- For a NW with N heterostructures, each of which in isolation 



p-5 



Matthew Shelley and Arash A. Mostofi 



Left Lead 


< Conductor 5 


Right Lead 






















III 1 1 1 1 


1 1 
















1 1 

A i B 


1 1 1 1 1 1 


1 1 

A B 





Fig. 4: Schematic illustration of MHNWs studied. Top: periodic conductor; middle: a typical 
random arrangement; bottom: Fibonacci chain pattern, with units A and B (see text). Si and Ge 
sections are in blue and magenta, respectively. 



gives a transmission Ti{E), we compute the transmission function as [35] 

r°'""(£;) = '^—^ , (3) 

^ ' N-{1- N)T/Tq ^ ' 

where T = '^iLiTl/N is the average transmission of the isolated heterostructures, and To 
is the transmission for the pristine NW. Having obtained T°^'^{E), k°^™ is calculated from 
Eq. ©. K™'^ and are used to estimate upper and lower bounds for zT. 

We consider MHNWs in the 111 and (211) growth directions, with diameters of 1.14 nm 
and 1.06 nm, respectively, and Gc heterostructure lengths of 3.80 nm and 1.34 nm, respec- 
tively. For each of these two growth directions, three qualitatively distinct heterostructure 
distributions are considered (shown in Fig. 21): (i) random, (ii) periodic, and (iii) Fibonacci 
chairl^l, each with ten heterostructures along the length of the NW. Structural details are 
given in Table [TJ For MHNWs with a random distribution of heterostructures, T(e) and 
T{io) arc ensemble averaged over 250 independent configurations. 

Fig. [5] (top panel) shows the thermal conductances of the MHNWs detailed in Table [1] 
and compares them to the pristine values. A prominent feature is the large reduction of Ki 
due to heterostructuring, with 111 MHNWs displaying smaller values than (211). Using the 
coherent model (red solid bars/blue bold striped bars), ki is reduced by factors of between 
five and eight when compared to the pristine results and reduces as the disorder is increased 
(periodic to Fibonacci to random patterning). The ohmic model (red shading/fine blue 
stripes), results in reductions of ki by factors of ~ 12 and 8.5 in 111 and (211) MHNWs, 



Fibonacci chain is an example of a ID quasicrystal 1361 : it displays local translational symmetries, 
yet remains aperiodic in toto. Exceptionally low k; values have been reported experimentally for 3D qua- 
sicrystals 1371 . thus the introduction a Fibonacci chain distribution of heterostructures could be a systematic 
method to reduce ft;. The Fibonacci chain MHNWs arc designed such that the length ratio of structural 
units A and B that comprise them is as close as possible to the golden ratio (1 + \/5)/2. Those structural 
units each contain a Ge heterostructure between lengths of SiNW and the total chain is built with three 
iterations (n = 0, 1, 2) of the sequence: A„^i = AnBn, B„+i = A„, with Ag = A and Bq = B. 





111 


(211) 


Single heterostructure length (nm) 


3.80 


1.34 


Periodic 


Total MHNW length (nm) 


93.3 


50.9 




Total number of atoms 


7208 


3608 


Random 


Total MHNW length (nm) 


93.3 


49.6 




Total number of atoms 


7208 


3520 


Fibonacci 


Total MHNW length (nm) 


93.3 


49.6 




Total number of atoms 


7208 


3520 



Table 1: Structural details of the MHNWs in the 111 and (211) growth directions. The different 
MHNWs (periodic, Fibonacci and random) are built by placing ten single heterostructures along 
the length of the wire. 
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Fig. 5: Transport properties at 300 K of 111 (red, solid bars) and (211) (blue, striped bars) periodic, 
Fibonacci and random patterned MHNWs, and comparison to the pristine cases. Results from fully 
coherent phononic transport 111 ((211)) calculations a have solid shading (bold stripe), while results 
in the ohmic regime have a lighter shading (fine stripe). Top panel: the lattice (left) and electronic 
(right) contributions to the thermal conductance, respectively. Bottom panel: zTmax (left) and S^G 
(right) at the value of fi that maximizes zT for each system. 



respectively. We note that Kf^^ / and K°'™/Ke are found to be between 2.5 and 7 — only 
marginally smaller than those values we obtained for single heterostructure NWs. 

■z^max and S'^G are shown in Fig. [5] (bottom panel) H It is striking that in no case 
does S'^G increase due to heterostructuring, and the (211) direction performs best over 
the range of MHNWs considered, showing only small decreases (with respect to pristine) 
as the disorder increases from periodic to random. Together with the pronounced effect 
that increased disorder has on k™*^, we see that, in the coherent regime, random patterning 
results in values of the figure of merit as high as zT = 2. Conversely, 111 MHNWs display 
significant reductions in S'^G as the disorder increases, which tend to counteract similar 
decreases in nf^^ leaving zT^^x approximately constant at ~ 1.7. In both 111 and (211) 
MHNWs the calculated zT^ax increases if the thermal transport is assumed to be ohmic, 
since «;°'™ < nf^. In this regime, nf^™^ is invariant with respect to the distribution of 
heterostructures, therefore, zTmax in this regime will follow the behaviour seen in S'^G 
with a value of zT ~ 2.3 for the (211) direction (almost independent of the distribution of 
heterostructures), and up to zT = 3 in the 111 direction, with a periodic arrangement of 
heterostructures . 

The effect of variability of heterostructure length. — The MHNWs that we have 
discussed thus far consist of multiple instances of identical heterostructures. Experimental 
synthesis techniques do not have this level of atomic precision, therefore, we have investi- 
gated the effect of introducing some variability of the lengths of the heterostructures that 
comprise the MHNW. In particular, we compare periodic patterned MHNWs with identical 
heterostructures with 'near-periodic' MHNWs that are comprised of heterostructures whose 
lengths are Gaussian distributed about a mean length that is given by the heterostructure 

^The electronic properties are always calculated within a fully coherent model but, depending on whether 
Kf'^™ or K^""^ is used, the value of /i at which ^Tmax occurs changes sUghtly and, hence, the value of S^G, 
which is dependent on fi. 
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Fig. 6: Comparison of S^G (left panel), k.i (middle panel) and zTmax (right panel) at 300 K in 
periodic (green triangles) and near-periodic (red squares) 111 MHNWs as a function of the number 
of Ge heterostructures (hets.). Results using the coherent (coh.) and ohmic (ohm.) phonon 
transport models are shown with solid and dashed lines respectively. 



length used in the 'true' periodic case, with a standard deviation a that corresponds to 
approximately 1/3 (1/2) a unit cell in the 111 ((211)) direction. For the near-periodic MH- 
NWs, the transmission functions are ensemble averaged over 250 independent configurations 
of the disorder to model a 'typical' MHNW of this type. 

We consider MHNWs with up to twelve heterostructures along their length. Fig.lHlshows 
a comparison between near-periodic (red squares) and periodic (green triangles) MHNWs in 
the 111 growth direction. The resulting values for zTmax are shown in the right panel. In the 
ohmic phonon transport regime (dashed lines), the near-periodic system displays a dramatic 
reduction in zT^ax as compared to the periodic case, which arises from the sharp reduction 
that is found in S^G (left panel), combined with the fact that k/ does not decrease very much 
(middle panel). When considering the coherent regime (solid lines), the reduction in S'^G 
for the near-periodic MHNW is also large (left panel) , as compared to the periodic MHNW, 
but associated decreases in m are also observed (middle panel) so that the resultant drop in 
•z^max (right panel), as compared to the periodic case, is much less pronounced. For the near- 
periodic MHNWs, both phonon transport regimes display a maximum in zTmax with respect 
to heterostructure length after the introduction of approximately four heterostructures. 

We note that the significant decreases in S'^G due to the variability in the heterostructure 
length is consistent with our earlier conclusion that increased disorder tends to reduce the 
power factor, as was seen when comparing periodic, Fibonacci and random distributions. 
We also note that in the ohmic phonon transport regime, there is little difference in ki 
between periodic and near-periodic MHNWs, which follows the earlier observation that 
there is almost no dependence of ki on heterostructure length in single heterostructure NWs 
(Fig. 131 top panel). 

Finally, comparing near-periodic MHNWs in the 111 and (211) growth directions, we 
find that S'2G'<"i> < 52G<2"), which may have been expected from the stronger dependence 
of S^Gniax on heterostructure length in 111 SiNWs with a single heterostructure (Fig. [5]). 
We also find that k^^^^ < which also could have been predicted from the trends 

observed for single heterostructure SiNWs (Fig. [3]). However, the delicate balance between 
S^G and ki make it difficult to use calculations on SiNWs with a single heterostructure to 
predict trends in zTmax for our near-periodic MHNWs, highlighting the need for accurate 
first-principles approaches. Among the near-periodic MHNWs studied, the 111 direction 
with four heterostructures performed best, with zT ~ 1.5 — 1.6. 

In conclusion, we have performed first-principles calculations on thin, p-type (110), 111 
and (211) SiNWs with Ge heterostructures. In all cases studied, a decrease of thermoelectric 
power factor S^G is observed when a heterostructure is introduced, and any increase in the 
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figure of merit zT is due to a corresponding reduction in k; . We have built model Hamilto- 
nians for MHNWs with over 8400 atoms while retaining first-principles accuracy. A similar 
method was applied to the dynamical matrices of MHNWs to obtain the thermal conduc- 
tance Ki for such structures. In such MHNWs we again find that S^G is always reduced 
and that increases in zT are driven predominantly by significant decreases in kj. We find 
values as high as zT = 3 in 111 MHNWs with periodic arrangements of Ge heterostructures. 
The intricate balance between S^G and kj, however, makes zT strongly dependent on the 
details of the system at the atomic level: in structures that model the kind of disorder that 
may be present in realistic MHNWs, more modest values of zT = 1.6 are found, which is 
still a factor of four greater than the pristine SiNW case. Our calculations suggest axial 
heterostructuring to be a promising route to high-zT nanowire thermoelectrics. 

* * * 
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